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I. INTRODUCTION. 

Cuprates in the metallic doping range have a large 
Fermi surface which is populated by short lived quasipar- 
ticles with lifetime strongly dependent on temperature. 
It has been proposed^ that short lifetime and strong tem- 
perature dependence are a result of coupling to fluctua- 
tions of a local charge state, which is independent of but 
coexists with the Fermi surface quasiparticles; to account 
for main experimental features the dynamics of the fluc- 
tuations of the charge state is to be local in space and 
non-local in time, Q(rt; r't') oc 5(r — r') x [t — i'] -1 ; this 
form of correlation function (dynamic mode) has recently 
been confirmed by microscopic calculation, ImQ(u) oc 
thw/2T; it obeys the uj/T scaling, typical of the dynam- 
ics of fluctuations near the quantum critical point.— On 
the microscopic level, the local charge state is a pattern of 
inter-atomic currents within copper-oxide unit cell 3 - (loop 
currents) which has a point group symmetry equivalent 
to a polar vector parallel to the copper-oxide plane, a 
toroidal moment. Once non-zero, the only dynamics 
possible for a local toroidal moment is to rotate its orien- 
tation. The part of the Hamiltonian which describes cou- 
pling between the Fermi surface quasiparticles and the lo- 
cal charge state has a form H C oupiing oc f/(r)[V x j(r)}. 
Here operator U(r) acts on the charge state; it generates 
rotation of the local toroidal moment around the axis 
perpendicular to copper-oxide plane^. The Fermi surface 
operator [V x j(r)] (where j(r) is quasiparticle current) 
has a symmetry of angular momentum perpendicular to 
the copper-oxide plane. Using a language of group the- 
ory, the dynamic mode QAuS) couples to Fermi surface 
quasiparticles via an operator which transforms under 
A-2 g representation of the lattice point group, the same 
representation to which z-axis component of the angular 
momenta belongs. 

Raman scattering directly probes charge dynamics 
in the system; the Raman cross-section da oc [1 + 
riB (w)]Im^„(u) is proportional to the Raman response 
function, Xa(u) = Xa(u,q = 0) = ((p a p a )) u where p a is 
charge density operator in the lattice symmetry channel 
a. In cuprates (which typically have tetragonal lattice) 
the cartesian product of the polarization vectors of in- 
coming and outgoing photons breaks into irreducible rep- 
resentations a = Ai g , A2 g , Big, Big, Eg, see, e.g., Ref. |6|. 
These are symmetry channels for which experimental 
data for Raman scattering in cuprates is specified. Nat- 
urally, Raman signal in Ai g symmetry channel would 



probe directly dynamics of the local charge state and re- 
veal the uj/T scaling, x( CJ )a=A 2g K thcj/2T. It is known 
that Raman spectra with similar behavior is observed 7 - 
in symmetry channels, other than Ai g . In this paper we 
demonstrate that Raman response in symmetry channels 
such as Big, Big originates on the Fermi surface; we show 
that particle- hole response exhibits uj/T scaling if quasi- 
particle coupling to the dynamic mode is allowed. The 
experimental situation will be addressed in a separate 
publication. 

The response functions are analyzed in terms of dy- 
namics of excitation states; in metals near the Fermi sur- 
face these are particle-hole pairs. In conventional metals 
the dynamic process is particle-hole rescattering medi- 
ated by screened coulomb interaction^. In cuprates in- 
teraction with the dynamic fluctuation has two major 
effects on the quasiparticle physics. First, charge fluctu- 
ations serve as a time-dependent environment which is 
capable of exchanging energy with the quasiparticle; as a 
result, quasiparticle lifetime becomes short and temper- 
ature dependent. Second, dynamic fluctuation can medi- 
ate particle-hole rescattering; unlike conventional metals 
where the particle-hole rescattering event is local in time, 
St ~ h/ep, in cuprates rescattering event can be highly 
non-local in time. 

Similar situation arises in many other strongly corre- 
lated metals^—. In this paper we focus on cuprates in 
the "strange metal" phase d 13 ' 14 , the region of phase dia- 
gram in the metallic doping range at temperatures above 
superconducting and pscudogap transitions, though we 
expect that our discussion is applicable more generally. 



II. RESPONSE FUNCTIONS ON THE FERMI 
SURFACE IN THE METALLIC CUPRATES. 

The dynamic mode in cuprates in the metallic dop- 
ing range is specified by its (retarded) correlation 
functioni^ 



4.9. = X6Q(u,q) R = {{U{rt)U\r't')% 



lmQ(uj,q) R = 



th. 2^ ' 



0. 



LJ > A 



(1) 



where U(rt) operates on the local charge state; it gener- 
ates rotation of the local toroidal moment (loop current 
state), see Ref. H for details. The cutoff A is deter- 
mined from ARPES data! 7 and is found to be ~ QAeV. 
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The spectral weight Xq can be estimated from the results 
of polarized neutron experiments^ Coupling of the dy- 
namic fluctuation to the Fermi surface quasiparticles is 
specified in terms of the vertex 




JoUk-k' 



(fc, fc')cj,c fc / 



(2) 



where 70 is the coupling constant; operators Ck are 
fcrmion quasiparticle operators on the Fermi surface. 
Function 4>(k, fc') is a normalized angular harmonics on 
the lattice which transforms under Ai g irreducible repre- 
sentation of the D^h group. From the space-group sym- 
metry perspective it behaves as ik x fc'; see Ref. \3 for 
more details and derivation of Eq. ([2]). The elementary 
particle-hole rescattering event mediated by the exchange 
of the dynamic mode is specified with 

' i ) = g(e - e'; fc, fc') = ToXovQ(e - s')\cf>(k, k')\ 2 , 



-< <- 



k' 



(3) 



where v is density of states on the Fermi surface (per one 
spin species) . We define lattice angular harmonics F a (fc) 
via the decomposition 

g{e-e l ;k,k')=Y,9c t Q{e-e')F a {k)F a {k l ) (4) 

The coupling constant 70, the spectral weight Xq and 
the Fermi surface density of states v are absorbed in the 
definition of dimensionless coupling constants g a and do 
not appear explicitly in the calculation and results. In 
cuprates the density of states on the Fermi surface is 
non-uniform. The momentum integration is performed 
with 



d 2 p' 



v / d£ 



2tt 



(5) 



where function v(9) has the symmetry of the lattice, 
i.e., it belongs to a scalar representation A\ g . ; in 
the following we use notation 0(9) = v(9)/vq where 
Uq = J (d9 /2tt)v(6) is average density of states. The an- 
gular harmonics F a (k) are normalized with : 



dO 
2^ 



0(0)F a (6)Fp(6) = 5 a0 



(6) 



Analysis of ARPES data^i indicates that the coupling 
constant in the density, a = A\ g , channel (which controls 
the inelastic part of quasiparticle selfenergy) is of order 
of unity; we will use notation g$ = gA lg ■ On symme- 
try grounds, within the point group D^h of the tetrag- 
onal lattice the angular decomposition of \(j)(k,k')\ 2 in 
Eq. ([3]) contains only three irreducible representations 19 ; 
these are symmetry channels B\ g , £?2 g and A\ g . For 



spherical Fermi surface gB ls /gA lg = 9B 2g /gA lg = -1/2; 
for Fermi surface in cuprates these dimensionless param- 
eters can be determined from the experiment. We note 
that decoupling of the energy and momentum depen- 
dences in g(e — e'\ fc, fc') is a consequence of the fact that 
in cuprates the correlation function of the dynamic fluc- 
tuation Q(u>,q) is local in spacei^ 

In conventional metals the low-lying excited states are 
represented by a collection of particle-hole pairs on top 
of the Fermi surface. Response functions are analyzed in 
terms of particlc-hole rescattering which in mathemati- 
cal description is represented by particle-hole ladder dia- 
gram. Free propagation of particlc-holc pair, represented 
by the particle-hole section, is interrupted by a rescatter- 
ing events mediated by Coulomb interaction and repre- 
sented by Landau parameters. An observation which is 
a basis for the whole picture is that particle-hole section 
is singular near the Fermi surface; the form of singular- 
ity is determined by continuity equation which follows 
from particle number conservation property^ By iterat- 
ing particle-hole sections the particle-hole ladder collects 
(and exhaust) all singular contributions. This way the 
mathematical formalism does not make use of any small 
parameter; it only relies on the long quasiparticle lifetime 
which is necessary for developing a singular particle-hole 
section near the Fermi surface. In conventional metal 
the only channel for relaxation of quasiparticle energy 
is by emitting particle-hole pairs which are rare close to 
the Fermi surface because of the available phase space 
is small. In the metallic state of cuprates there is a new 
channel for quasiparticle energy relaxation. Emission and 
absorption of dynamic fluctuation is not confined by any 
phase space restrictions; in result there is a large inelastic 
quasiparticle relaxation rate near the Fermi surface and, 
consequently, short inelastic lifetime^ ' 13 i 14 

We base our discussion of the response functions in 
the metallic state in cuprates on the assumption that 
even though the particle-hole section has weaker singu- 
larity (branch cut, not a pole), it is still singular enough 
to justify the use of particle-hole ladder near the Fermi 
surface. Physically, this corresponds to analysis in terms 
of particle-hole rescattering. The technical side of the 
problem is more complex here compared to conventional 
metal. Indeed, particlc-holc rescattering mediated by 
exchange of dynamic fluctuation is non-local in time, 
Q(tt') ocl/t — t'. In the mathematical analysis we make 
sure that particle number conservation is obeyed. The 
calculation is controlled by maintaining the balance be- 
tween an inelastic quasiparticle relaxation rate ("outs") 
and particle-hole rescattering ( "ins" )&22r^2. In that sense 
the present analysis follows the spirit of the conventional 
Fermi liquid theory. 

Particle-hole rescattering physics near the Fermi sur- 
face determines only the dynamic part of the response 
functions . 15 i 20 This is enough for the purpose of calcu- 
lating the Raman response which requires calculation of 
ImXa(w); the static part of Xa( w ) i s reai - Complete 
particle-hole rescattering process consists of an arbitrary 
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sequence of two kinds of rescattering events : those which 
are local in time (represented by Landau parameter B a ) 
and those which are nondocal in time. In the particle- 
hole ladder we identify segments Z a (w) (singular seg- 
ment) which describe sequence of particle-hole rescatter- 
ing mediated by dynamic mode only. The correlation 
function x Q (w) has a structure of a geometric series in 
which there is an arbitrary number of singular segments 
interrupted by the Landau amplitude B a 



1 - B a ZS 



(7) 



Here cf> a is static Raman vertex (which also includes static 
Fermi liquid renormalizations). Landau amplitudes are 
defined via the angular decomposition of the non-singular 
scattering amplitudes B(k,k') = £ B a F a (0) B F a (0') B ; 
here F a {9) B are normalized angular harmonics. Index B 
indicates that these are different from angular harmonics 
F a {&); in the same representation their overlap is not 
equal one, f(d0/27r)P(0)F a (0)Fp(0) B = f a S af} , though 
they are orthogonal when in different representations. 

Analysis of particle-hole rescattering within the singu- 
lar segment leads to the following expression 

z Z = ~j «fei[thi±^-th|] 

x J^F a (0) B Y«(e) e F a (e) B (8) 
where the function Y(e) is determined from the equation 



uY„(e)<§ = 1 



v{0')d0' 
2^ 



g a F a {0)F a {6')YZ{e') e > - g o F o {0)F o {0')Y«{e)e 

"(9) 

where kernel Q% (e' — e) is defined in Eq. (|52"j) ; see next 
section for details. The singular segment is interrupted 
by non-singular rescattering event at both ends, hence 
the angular harmonics F a (0) B in Eq. ©. In the density 
channel, a = A\ g , the solution of Eq. (J5J is Y u {e) = 1/uj 

and consequently Z^ 9 = — 1. This ensures particle 
conservation; see detailed discussion in Section IIII CI 
Eq. (|60p has a form reminiscent of the kinetic equation, 
the integral on the right hand side plays a role of the 
collision integral. It is shown in the next section that 
the term — goY"(e) originates in the inelastic quasiparti- 
cle scattering rate ("outs") and the term g a Y"(e) origi- 
nates in the process of particle- hole rescattering ("ins"). 
Eq. replaces kinetic equation in the metallic state 
of cuprates; we have to restore the g-dependent term in 
the left-hand side, account for elastic disorder and local 
in time rescattering events, and replace the trigonomet- 
ric factors in the definition of Qff by the fermion and 
boson distribution functions. In this form Y u>q [e) now 



depends on q and is suitable for the analysis of transport 
properties; these however require separate investigation. 
We note that the structure of Eqs. (|6QI6ip is indepen- 
dent of the particular form of correlation function of the 
dynamic mode. They are valid for a description of the 
particle-hole rescattering near the Fermi surface in the 
presence of dynamic mode of arbitrary form if validity of 
particle-hole approximation has been established. 

Analysis of Eq. © in Raman channels leads to the 
result 



(l + B ot )uj + 2(g -g ot )^{w) 1 



(10) 



The function ^(w) has the same analytic properties as 
quasiparticle selfenergy; it has a form 



* Q (w) T = 



2ds 



* q (cj-s) t=0 

Im* Q (w) T =o =|w|, 

Retf Q H T =o =wln^ ^ 
M 



where p « 0.59. The cutoff A a cannot be determined 
within the present calculation because it is not deter- 
mined by the physics in the vicinity of the Fermi surface; 
it is reasonable to expect that A a is of order of the high 
energy cutoff in the fluctuation mode, Ao- Note the fac- 
tor g a — go- It is a result of partial cancellation between 
the dynamic singularity in the selfenergy (oc go) and the 
singularity introduced by non-local in time particle-hole 
rescattering, both are due to the coupling of the Fermi 
surface quasiparticles with the dynamic mode. The ana- 
lyticity of Xa( w ) m the upper half of the complex plane is 
guaranteed by the fact that g a — g > 0. The imaginary 
part of response function in Raman channels a ^= A\ g 
calculated from Eq. (|T0|) obeys oc u/T scaling at oj < T; 
at large energy, uj >T, it levels to a constant. 

In cuprates the dynamic mode can only induce 
particle-hole rescattering in "tensor" symmetry channels, 
Ai g , Big, B2 g \ it does not participate in the particle-hole 
rescattering in the "vector" channels, such as E g and E u . 
We now briefly discuss consequences of this fact. The 
current vertex (charge, heat, etc.) is spatial vector and 
it belongs to E u irreducible representation. Therefore 
in the direct calculation of the diagram for conductivity 
tensor the dynamic mode enters only via selfenergy; no 
particle-hole rescattering mediated by the mode is pos- 
sible. This justifies single-bubble calculation of conduc- 
tivity. The case of long-range disorder is captured by 
clastic disorder ladder in the transport diagram; this is 
equivalent to the discussion in terms of conventional ki- 
netic equation^. An example of the transport property 
where particle-hole rescattering via dynamic fluctuation 
is important is Fermi surface contribution to the solid 
viscosity^ which enters in the acoustic sound attenu- 
ation. Indeed, the stress tensor in the tetragonal lat- 
tice breaks into Ai g , Bi g , Big, E g irreducible pieces. The 
solid viscosity is determined by the dynamic stress-tensor 
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response function. We note that, unlike Raman corre- 
lation functions for which qvp <C w, in acoustic sound 
the relation between frequency and momenta is reverse, 
qvp 3> cj; this requires analysis of q-dependent kinetic 
equation; see discussed after Eq. (|9|). 



III. PARTICLE-HOLE RESCATTERING 
MEDIATED BY DYNAMIC FLUCTUATION 

A. quasiparticle selfenergy 

In the remainder of this paper we indicate sequence of 
mathematical steps necessary to arrive to Eqs. (|9I8I10[) . 
We use temperature technique which allows direct access 
to analytic structure of physical objects at hand; for de- 
tails on analytic continuation in temperature diagrams 
see, e.g., Ref. [2ft To make discussion more readable we 
assume the Fermi surface is spherical; angular harmonics 
can be easily restored as is discussed above. 

Emission and absorption of the dynamic mode by the 
quasiparticles on the Fermi surface results in an inelastic 
relaxation rate which is analyzed in terms of the selfen- 
ergy : 

iuj;p— p' 



E(«e„) = .? , f ) \) 



d 2 p' 



9oT^V] I —j^—Tn G(ie n - iu n )Q(iuj n ) 



Making analytic continuation we obtain 



(12) 



f \ l£n>0 



ie ie—iu) %e 



d 2 p' 



coth 



R—A 



^o(2tt) 2 

R 



2T R 



fcr — b 



th- 



2T 



are G R G R or G A G A pairs. The energy integration is not 
confined to the Fermi surface here; naive momentum in- 
tegration for the pair of Green's functions gives zero since 
in the complex plane of £ both poles are on the same side 
and expression falls fast enough ~ l/£ 2 at large £. How- 
ever, this contribution is not zero: as noted above the 
frequency integral is not confined to the Fermi surface in 
RR section; the combined dsd£/[(e — lo — £)(e — £)] inte- 
gral is logarithmically divergent at the upper limit and as 
usual have to be cut off at the Fermi energy scale. How- 
ever, we observe that this contribution does not lead to 
a singular behavior in e (or an imaginary part). We only 
collect the singular (non-analytic) at the Fermi surface 
terms. 

Momentum integration is trivial in Eq. (|14p owing to 
the fact that the mode is local in space; we use an identity 



/ 



d 2 p' 



ImG R (s,p) 



d£lmG R (s,0 



-VqTX 

(16) 



and obtain 



ImE*( £ ) = 90 J dw±[coth^ - th^] lmQ R (u>) . 

(17) 

The real part is calculated from Kramers-Kronig rela- 
tions 



Re/H = H.T. {Im/H} = dio'^f^ 

We use Eq. (pj in Eq. ^TTJ) and obtain (at e <C A ) 
ImE(e) = - g - 



(18) 



'th[e/2!T 
ReE(e) = - g R.T 



th[e/2T] 



(19) 



To analyze this expression we introduce a representation 
which express the finite-temperature self-energy in terms 
of its form at zero temperature. To illustrate the idea we 
first calculate Q(lu)t- Introduce a "temperature smear- 



(13) ing" operator- 



26 



or 



^r, \ f dto f d 2 p' 



x [coth ±L - th^]lmG R (e - u,p')lmQ{u) (14) 



Starting with 



2ds 



ItciQ,(uj)t=q = — signw 



(21) 



Let us comment on the higher order terms in the selfen- 
ergy. 



and using an identity th(w/2T) = $ signo; fi„ =1 , we write 
the imaginary part of the selfenergy, Eq. ([Tj) , in the form 



) A f > 



(15) 



Making an analytic continuation in this diagram we find 
that all pairs of Green's functions, except the middle line, 



ImQ R (u) T = a ImQ fl (w) T=0 8 P=1 



(22) 



An important property of the operator 8 ■ • • 11 p is that it 
commutes with Hilbcrt transform if we restrict ourselves 
to functions which are analytic in the upper/lower half 
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of the complex plane. Indeed, smearing operation repre- 
sents function g(oj) as a linear combination of functions 
f(oj) shifted parallel to the real axis. If function f(oj) is 
analytic in the upper half of the complex plane, so is the 
linear combination in the right-hand side of Eq. (|20|) . As 
a consequence, both sides in Eq. ([20| satisfy Kramers- 
Kronig relation with branch cut along the real axis. We 
immediately obtain 



RcQHt = 8 ReQ(w) T=0 8, 



and therefore 



QMt = J ds (^^) C(w - s) T=0 



At T = one have 

ReQ(w) T =o = H.T. {RcQ(w) T=0 } 



2. A 
-ln—r 

7T \U>\ 



(23) 



(24) 



(25) 



Using smearing representation in Eq. (|24|) we write 
fdth-, 



ReQ R (uj) T = / ds 



l 2T 



V 2ds 



2. A 
-In, 1 



(26) 



Expression in Eq. (|26[) gives a very intuitive representa- 
tion for Q(oj)t at finite temperature. It makes it trans- 
parent that when oj 3> T the difference between Q(uj)t 
and <2(w)t=o is small. It also isolates the logarithmic di- 
vergence in the Hilbcrt transform which exists already 
at zero temperature. At small oj <C T the integral 
in the real part can be evaluated in closed form using 

fda(f%p)}n\8\ 
a finite limit as oj 
formula 



—7 — ln(4p/7r); real part tends to 
0. We can write an interpolating 



RcQ r {oj) 



it max{w, C 1 \ 



where C = ne^ /2 « 0.88 and 7 « 0.57 is Euler's 
gamma. 

We now obtain similar "smearing" representation for 
the quasiparticle selfenergy at finite temperature. At zero 
temperature the selfenergy, Eq. (fT?)) . is easily evaluated 



lmT, R (e) T =o = So 



The real part is 



/ dojImQ R (oj) T =o = -5 min(|e|, A ) 
Jo 

(28) 



ReE(e) T =o = H.T. {ImE(£) T =o} 
2 



-0o- 

7T 



" A 1, . e+A 1 e-A - 
e ln_ + -( £+ A )ln— — + -(e-A In— — 
el Ao 2 Aq . 

2 r . A 



- 50 - [e In ^+6] (29) 



Following the analogous discussion for the dynamic mode 
Q(w) we first find the smearing representation for the 
imaginary part. Using an identity 

2- ui-s = -In 2 cosh L — 30 

2as p L 2 J 



(is 



and comparing with the finite temperature expression for 
the self-energy, Eq. ([IT)]) , we write 



2T 

Im£ fl (£)«-<7oyln 



2 cosh 



pe 
2T 



(31) 



For any value of p this expression has correct behavior 
at e 3> T. The choice of p is determined by the form of 
selfenergy at small e < T. One choice is p = In 2 « 0.69 
which ensures correct value at ui = 0. The other choice 
is p = 2/3 ~ 0.66 which ensures correct value of the 
curvature at s = 0. Two values of p are close to each 
other and indeed for any value of p between these two the 
difference between two sides of Eq. (|5T|) is barely visible 
on the plot. Eq. pip provides a good interpolation in the 
whole range of e. We conclude 



IhiS(w)t = 8 ImS(cj)T=o 8 



and therefore 



E(w)g 



/dthf| 
I 2ds 



E(w — s) 



T=0 



(32) 



(33) 



with p w ln2. Note that p = 1 in Eq. (|M|) . 

Finally, we indicate how these results arc modified 
when density of states is non-uniform on the Fermi sur- 
face. Eq. (|14p is replaced with 



Im£( £ )f = -g F {9) 



2^ 



d 2 p' 
(2tt) s 



x [coth ^ - th^— %mG*(£ - W ,p')ImQ(a 



2T 



2T 



- 9o F (6)(F )ImE(e) R (34) 



where ImE(£) fi is given by Eq. p7|) . Here (i 7 !)} = 
/ d6(2ir)i)(6)Fo(6). For non-circular Fermi surface the 
constant go (Fq) plays the role of the dimcnsionless cou- 
pling constant in the density channel. The selfenergy 
acquires slow dependence on the position on the Fermi 
surface^. 



B. Calculation of singular segment Z". 

The singular segment Z™ describes multiple rescatter- 
ing of particle-hole pairs via an exchange of a dynamic 
mode 



Z a = 



(35) 



With our definition of the singular segment, Eq. ([7|. 
the triangles in this expression correspond to a factor 
of one. In Ref. [IH Zo(oj) has been expressed in terms of 
a fully rcnormalizcd particle-hole propagation amplitude 
r 7 / /7 (£i, £ 2 ; oj) for which (Bethe-Salpeter-like) equations 
have been derived. In constructing these equations the 
analytic structure of r / / // (£i, £2; oj) is essential and it 
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was discussed in detail. Here we take simpler approach 
which allows for a more transparent discussion. We note 
that all major results can also be obtained by continuing 
the line of analysis in Ref. HH 

We start by calculating each term in the series for 
Z° in Eq. ([55)1 . We introduce a function Z"(e) which 
is obtained by omitting the last section and the corre- 
sponding energy sum in each term in Eq. (|35j) : we define 

Z u (e) — X)n=o oo Z& (e) (to avoid excessive indexing we 
omit a index where it does not lead to confusion). Zero 
term is one 



zL 0) (s) = i 



The first term is 



Sl+UJ 

T 



(36) 



(37) 



The function Z u {e) has only three analytic pieces in the 
complex plane of £2; these are separated by branch cuts 
at I me 2 = and I me 2 = —u n , i.e., its analytic structure 
is equivalent to that of the particle hole section which 
describes free propagation of particle-hole excitation. To 
label the analytic pieces of Z^{e) we use notation which 
is commonly used for particle-hole section : 



A A — > Ime < — oj n 
RA -> - uj n < Ime < 
RR^ < Ime . 



(38) 



We have assumed that Imu; > 0, i.e., we are calculating 
retarded part of the response function. The main advan- 
tage of working with Z u [e] is that in the RA range of ie n 
it has only one analytic piece. Compare this with the full 
rescattering amplitude Y 11 / 1 analyzed in Ref. UH which 
has two analytic pieces in the RA range of frequencies. 

The analytic continuation for the RA analytic piece 
has a structure 




-th — \n 
2 2^JU 



1 



coth 



£1 - £2 



> , > 



< ' < 



ii-i 



(39) 



Here and in rest of discussion we omit temperature in the 
trigonometric factors where it does not lead to confusion. 
Note that the amplitude for the rescattering event (given 
by the vertical wiggly line), viewed as a function of e\, £2, 
has the same analytic structure as the fully renormalized 
amplitude r 7 / 7/ (ei, e 2 ; w); since only one rescattering 
event is involved here, the analytic pieces //// directly 
correspond to RJA analytic pieces of Q(w). The first 
and second terms in Eq. (|39[) arc generated by the branch 



cuts at Imei = and Imei = — uj n . The third term orig- 
inates in the dynamic mode which introduces a branch 
cut at Imei = Im£2. We rearrange Eq. (|3^|) in terms of 
Q+ = Q" + Qj = 2ReQ and Q~ = Q u - Q 1 = 2TmQ 



Zi 1] {e 2 )RA = 9a dei 



1 



- +2 coth - 



SII+I 



> , > 



< ' < 



(40) 



In this form the energy integration is restricted to the 
vicinity of the Fermi surface in both terms. The mo- 
mentum integration in the RA particle-hole section is 
done by deforming contour in the complex plane of mo- 
mentum represented by an energy variable £ defined via 
v£ = e\ p +8p[. We define a function S(u>)ra 

/j2 
^JL.G* +cl (e + u)G*{e) 



d6 



2-7T uj — qvp cos 6 — [S fl (e 



X A (e)} 



(41) 



In the complex plane of £ the poles of two Green's func- 
tions in the RA section are on different sides of the con- 
tour of integration, the distance of the pole from the 
real axis is proportional to the imaginary part of selfen- 
ergy with the coefficient 1 + dReE/dfc. Since in cupratcs 
the momentum dependence of selfenergy is very weak, 
when imaginary part of selfenergy is small, the pole in 
the complex plane of £ is close to the real axis. This is 
the Migdal theorem^ applied to cuprates. ^-integration 
extends to a large momenta, £ ~ ep, therefore for e,u> 
close to Fermi surface the momentum integration can be 
done by taking a residue at the pole on one of the sides 
of the real axis (the contour pinching). Though it is not 
obvious that in cuprates (where the Green's function has 
a singular selfenergy) analytic structure of the particle- 
hole section in the complex plane of £ is similar to the 
conventional case, the calculation path taken here is con- 
sistent because we make sure that the conservation laws 
are obeyed (see below). For non-spherical Fermi surface 
Eq. (|4Tj) is replaced with 



2 w 5(e; W ,< 7 ) Q = J ^JLF2(p)GK +q (e + u)G£(e) 
P(6)d6 Fl(9) 
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- qv F cos0 - [E(e + u)f - T,{e)f] 



(42) 



Analytic continuation in RR and AA analytic pieces 
of Zu (e) leads to an expression in which the energy inte- 
gration is no more restricted to the vicinity of the Fermi 
surface. The combined energy and momentum integra- 
tion is logarithmically divergent, compare this discussion 
with the one following Eq. (fT5|) . However there is no sig- 
nificant contribution to the non-analytic behavior near 
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the Fermi surface since both poles in the complex plane 
of £ are on the same side of the contour of integration; we 
only consider RA analytic piece in the ladder (and omit 
RA subscript). 

To keep discussion simple and to illustrate the main 
point we first restrict ourselves to small frequencies, 
uj <C T and omit all terms that vanish in this limit; 
this means we omit the real part of the selfcnergy and 
the first term in Eq. (|40|) . It will be trivial to restore 
complete equations later on, see discussion leading to 
Eqs. (JgDJ and (JBTJ. We write 



Z£\e 2 )=g a f 



-th-^-f coth : 



> , > 
11-1 



<fei^[coth £l - " 2 



th £ j-]S(e 1 )Q-(e 1 -e 2 ) 



r J[l-th^th^] (43) 



uj + 2ig e coth ^ 2 



where 



To obtain we have to add the last section and perform 
the frequency sum. Here the analytic continuation in the 
frequency sum leads to a standard trigonometric factor: 
in the RA section there are no additional branch cuts. 
We obtain 



Z, 



deZJe) --[tbJ^-^ - th-1 

" V w + 2i 9o ecothf 2 L 2 2 J 



dsZ 0J (e) 



1 



uj + 2igo£ coth ■ 



uj- 



dthf 
2de 



(49) 



the minus sign is a result of our definition of Z u (e) which 
starts at Z^ (e) = 1; to check the sign consider g a = in 
which case calculation of Z^ reduces to a single bubble 
diagram with the result — 1 . In the second line of Eq. (|4"9")l 
we have used 



th- 



+ uj e 

th- 

2 2 



-¥ UJ- 



cfthf 
2de 



(50) 



C _ (w)=-2ith— 



and 



S{s) 



1 



uj + 2igoS coth ■; 



Now consider the next term in the series 



? 



(44) 



(45) 



(46) 



To make analytic continuation for the frequency sum in 
the second section we have to consider the analytic struc- 
ture to the left, Z^ 1 \e), and to the right, Q~(e2 —£3). It 
has been noted that Z u (e) has only one analytic piece in 
the RA range of Matsubara frequencies. We only have to 
consider a branch cut at Im£2 = Im£3 which originates in 
the second rescattering event and therefore the analytic 
continuation for the frequency sum is done the same way 
as in the case of Z^ x \e). Keeping only terms that do not 
vanish in the limit uj — > 0, we write 



^ (2) (£ 3 ) 



1 9 9 9 1 



uj + 2igoei coth ^ 2 
2ig a 



uj + 2igoS2 coth 4f- 2 



The rest of the terms Zu(e) are constructed the same 
way. The form of each term in the series Z^\e) indicate 
that the sum Z^(e) satisfies an equation 

Z u (e) = 1+ 

fds'Z u (s') ^ _I[i_th^-^tl4] (48) 

J uj + 2igoe' coth ^2 2 2 



^l-th^^th^] (47) 



C. particle number conservation. 

Particle conservation requires that the density correla- 
tion function vanishes identically for q = 0; it is therefore 
necessary that in this limit the dynamic part of the den- 
sity correlation function approach real, oj-independent 
value. With our definition of the singular segment we 
have to check that Z^~ Alg = — 1 if calculation is made 
in the density channel, A\ g ; see Refs. [H and I20I for de- 
tails. We start with Eqs. (|4"5)l and (|4^1) which have been 
derived for uj <^ T . Eq. (|49| has an overall factor uj; in 
the limit uj — > the function Z^ will vanish unless Z u (e) 
is singular; we conclude that to ensure particle conser- 
vation function Z ul (e) has to behave as Z u {e) oc 1/oj 
at small omega. To see that this is indeed the case, 
we observe that at finite temperature the denominator 
of the particle-hole section S(u) is finite and the fre- 
quency sum in Eq. (|40|) has a contribution which does 
not vanish in the limit uj — > 0. Therefore each diagram 
■2T(i=o( e ) is finite (and positive) and at uj = the sum 



Z^2o( £ ) wm diverge. We expect that the 



Z w =o(e) 

effect of a small but finite uj is to regularize the sum with 
the result Z u (e) oc 



We see that the fact that Z^2o( £ ) i s finite in the limit 
of uj = is crucial for the consistency of our calculation. 
The trigonometric factor in the first term of Eq. (|4"0)) 
vanishes in the limit of uj — > while the second remains 
finite. Physically, this nonvanishing contribution origi- 
nates from the singularity in the dynamic mode, the same 
singularity which is responsible for the inelastic behavior 
of the quasiparticlc selfenergy. There is a parallel with 
the calculation of diffusion ladder for the elastic disor- 
der; in the disorder ladder the absence of a factor of uj in 
each particle-hole section is natural because there is no 
frequency sum in the intermediate sections. In the case 
of dynamic mode the finite value (in the limit uj — > 0) of 
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imaginary part of sclfcncrgy ( "outs" ) is compensated by 
the finite result of the frequency sum in the particle-hole 
section ("ins"). Similar to elastic disorder, their complete 
cancelation in the calculation of Z"~ Ala ensures particle 
conservation. Unlike disorder, it does not seem possible 
to calculate each term Z& (e) separately and follow the 
cancellation explicitly. 

We now observe that Eq. (J3HJ) has a closed solution in 
the density channel : 



Z u (e) 



UJ 



2ig e coth ■ 
2ig 



(51) 



To check, we substite this expression into the right hand 
side of Eq. 



Z u (e) = 1 + 



2ig 



th- 



-£,£'. 

th— 

2 2 



(52) 



and use the fact that the integral here has the same form 
as the one which appears in the calculation of selfcnergy, 
see Eq. ([17]) 



de' QT {e'-e) \ [coth f^-f-th-] 



2ig e coth ■ 



(53) 



We use Eq. (|53j) to rewrite it as 

w y°( e ) = 1 



de' 



2ig a Y«(e')-2tg Q Y«(e)\-ll 



th 



e' — e , e' 
2 2 J 



(58) 



In this form it is obvious that in the density channel, 
9a — go, the solution of this equation is independent of e 
and therefore has a form 



Y u (e) 



1 



(59) 



When substituted in Eq. (|56|) it gives Z^ = — 1. 

To obtain an equation for Y(e) which is valid for arbi- 
trary lu/T we start with Eq. (|4T)|) and repeate the same 
sequence of mathematical steps. The result is 



coYUe) a = 1 



de' 



g a Y«(e') - g YZ(e) Q%(e' - e) (60) 



which together with 



— -th-] (61) 



We have 



Z u (e) = l + 



2igoe coth ■ 



(54) 



which coincides with Eq. (|5Tj) . Substituting Eq. (|51|) into 
Eq. (J49J) we obtain 



7 a=A lg _ 



= -1 



(55) 



which proves particle conservation. 

In the rest of this section we extend the above analysis 
for arbitrary uj/T and restore the real part of the selfen- 
ergy (and the first term in Eq. (|40|)). First, we rewrite 
Eq. P8")l in the form which make the cancelation more 
transparent. Introduce a new variable Y°{e) which cor- 
responds to adding a particle- hole section, S(e) to each 
diagram in zj^'ie) 



determines the singular segment for arbitrary u and T. 
Here we have defined a function 



, £ ' - £ 1 r , S' 

coth— — --[th- 



th-] 



Q + {e' 



th- 



th- 



(62) 



In the density channel the solution of Eq. is inde- 
pendent of the temperature and coincides with Eq. (|59|) . 
This concludes the proof of a particle conservation prop- 
erty within the framework of our mathematical discus- 
sion. Eqs. (|60| and (|61[) are the main technical result of 
our discussion; for the purpose of the analysis of the re- 
sponse functions on the Fermi surface they replace ladder 
equations derived in Ref. [TH 



y°(e) =ZZ(e) 



uj + 2ig e coth | 
dtM 



deY«(s) 



u>- 



2de 



Y u] (e) satisfies an equation 



[w + 2i fl0 ecoth-]y°(e) = l 



(56) 



£ — £ 



+ 2ig a / de'YZ{e')-\l - th— -th-] (57) 



D. Calculation of Z" in Raman channels, g a 7^ po- 
ll! Raman channels, g a ^ go, the solution of Eq. (|60[) 
is a function of e and it appears that closed solution of 
the equation is not possible. Here we obtain approxi- 
mate solution which is consistent with the the analytic 
structure of the singular segment in Raman channels. We 
solve Eq. (|60|) for u) <C T and u ^> T and give and inter- 
polating formula which bridges the two cases. Analysis 
of experimental data will require more rigorous analysis 
of Eq. First consider the case of large frequencies, 
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u> ^ T. We omit the Q + part in the expression for Q K ; 
Eqs. (jnnj) and (EB) take a form 

max{— u,0} 

uY° (e) = 1 + 2^ J de' (g a Y£ (s') - g Y^ (e)) 

min{— cj,0} 
max{— o;,0} 

Z« = -signw /" del^(e) (63) 

min{— u,0} 

In the first equation we have used the fact that —lu < 
e < 0, see Eq. (|6"Tj) . The solution Y(e) is independent of 
e ; the function Y" satisfies an equation 



We find 



uY° = l + 2z(g a -g )\uj\Y° 



Z°[ = 



lu + 2i(g - g a )\u\ 



(64) 



(65) 



With this expression in mind we can write singular seg- 
ment in the form 



Z a , 



us + 2 \g - g a 



(66) 



where ^S(lu)t is analytic in the upper half of complex 
w-plane (i.e., it is a retarded function, same as Z(uu)). 
Comparing with Eq. (|65|) (and using results in the rest 
of this section) we write 



Im*(w)a 



CqT, 



lu >T 



(67) 



where Co ~ 2.3. Note that analyticity of Z u (lu) in the 
upper half of the complex plane is guaranteed by the fact 
that go — ga > 0. Having in mind our experience with 
the calculation of the selfenergy we write an interpolating 
expression 



, . r /dthffA 
Im*(o;)r = / ds I j!J jlm^o; - s)q 



(68) 



Where constant p is fixed by the behavior of $(lu) at 
lu < T; we find p « 21n2/2.3 0.59. The real part 
of ^(w) is calculated from the first line in Eq. (|67|) via 
Kramers-Kronig integral. Wc find 



Im$(w)r = o = Mj Rc v E'(w)t=o = win 



A 



*(w)<] 



I 2ds 



- s)t=o 



(69) 



The cutoff has to be of order of the upper cutoff Ao 
in the dynamic mode, see Eq. ([T]). 

In the rest of this Section we analyze Eqs. (|60|) and ([6"Tj) 
in the limit lu <C T which determines constant cq in 
Eq. ([B7|. Eqs. (JBDJ* and (JMH take a form 



w y°(e) = 1 

+ 2zJ de'(g a Y:(e')-goYZ{e))\ 
f dth§ 

Introducing new variables 



1 - th th— 

2 2 



(70) 



77 = th— , 7]' = th — . 

' 2' ' 2 
we rewrite Eq. (|70[) in the form 

<W =1 + 2% j dr) ± i- 



(71) 



1 ?^(?/) 



(72) 



Multiplying the first equation by J di^/2 we write 



1 d?7 



w + 2?(ff ~ 9a) 



2atanh?7 



V 



Y«( V ) = 1 (73) 



where we have used J", dr?'/(l — r/7/') = (2atanh?7)/?7. 
The solution is 



1 + c(r?) 



cj + 2i(.(/o - 9a) 



2 at an hn 



(74) 



where 0(77) is a function which has zero mean, J c(t/) = 
0; it can be determined by substituting Eq. (|74[) into 
Eq. (|72p; for our approximation we take 0(77) = 0. Sub- 
stituting this into Eq. (|T2")) and taking the limit lu — > 
we find Z a {ui) = —Lu/[—2i(g a — <7o)cor] where 1/co = 
/_, 77^77/ (4atanh?7) w 0.43 and therefore Co ~ 2.3. 
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